The Structure of the Local Interstellar Medium. III. Temperature 

and Turbulence 1 

Seth Redfield 2 , & Jeffrey L. Linsky 
JILA, University of Colorado and NIST, Boulder, CO 80309-0440 

ABSTRACT 

c ■ 

We present 50 individual measurements of the gas temperature and turbulent 
velocity in the local interstellar medium (LISM) within 100 pc. By comparing 

CN| \ the absorption line widths of many ions with different atomic masses, we can 

satisfactorily discriminate between the two dominant broadening mechanisms, 

>■ ' thermal broadening, and macroscopic nonthermal, or turbulent, broadening. We 

^ ■ find that the successful use of this technique requires a measurement of a light 

ion, such as D I, and an ion at least as heavy as Mg II. However, observations 

Q ' of more lines provides an important consistency check and can also improve 

the precision and accuracy of the measurement. Temperature and turbulent 
velocity measurements are vital to understanding the physical properties of the 
Q_i' gas in our local environment, and can provide insight into the three-dimensional 

morphological structure of the LISM. The weighted mean gas temperature in 
the LISM warm clouds is 6680 K and the dispersion about the mean is 1490 K. 
The weighted mean turbulent velocity is 2.24 km s" 1 and the dispersion about 
the mean is 1.03 km s -1 . The ratio of the mean thermal pressure to the mean 

^\ • turbulent pressure is Fr/Pg ~ 26. Turbulent pressure in LISM clouds cannot 

explain the difference in the apparent pressure imbalance between warm LISM 
clouds and the surrounding hot gas of the Local Bubble. Pressure equilibrium 
among the warm clouds may be the source of a moderately negative correlation 
between temperature and turbulent velocity in these clouds. However, significant 
variations in temperature and turbulent velocity are observed. The turbulent 
motions in the warm partially ionized clouds of the LISM are definitely subsonic, 
and the weighted mean turbulent Mach number for clouds in the LISM is 0.19 
with a dispersion of 0.11. These measurements provide important constraints on 
models of the evolution and origin of warm partially ionized clouds in our local 
environment. 
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1. Introduction 

Understanding the physical state of the interstellar medium (ISM) is a fundamental 
area of research in astrophysics. Theoretical studies of the phases of the ISM have produced 
classic papers in the astrophysical literature (e.g., Field, Goldsmith, & Habing 1969; McKee 
& Ostriker 1977), the details of which are still being discussed, analyzed, and refined to 
this day (Cox 1995; Wolfire et al. 1995a,b; Heiles 2001a). These classical theoretical models 
predict the stable coexistence in pressure equilibrium of various phases of interstellar matter: 
the cold neutral medium (CNM), the warm neutral medium (WNM), the warm ionized 
medium (WIM), and the hot ionized medium (HIM). Detailed observations are required to 
test the theory and provide direction for future modifications. 

The heliosphere is surrounded by warm partially ionized material (Redfield & Linsky 
2000) in the local interstellar medium (LISM), which has properties intermediate between 
those of the WNM and WIM theoretical models. The LISM is embedded in a larger cavity, 
the Local Bubble, presumably filled with hot rarefied gas. Very little cold material, or 
CNM, resides in the LISM, as is shown by the lack of Na I absorption towards nearby 
targets (Lallement et al. 2003; Sfeir et al. 1999; Welty, Morton, & Hobbs 1996; Welsh et al. 
1994). However, observations of this cold material are made from low ionization absorption 
line measurements (Jenkins & Tripp 2001), and from absorption and emission in the 21cm 
hyperfine structure line of H I (Heiles 2001b). The HIM that presumably pervades the Local 
Bubble, was identified by the presence of soft X-ray emission by Snowden et al. (1998). The 
Cosmic Hot Interstellar Plasma Spectrometer (CHIPS) instrument, which was launched on 
2003 January 12, should provide important new observations of this rarefied gas (Hurwitz 
& Sholl 1999). 

Observations of the warm partially ionized gas are vital to unlocking the physical struc- 
ture of the ISM. Distant WNM environments can be probed by analysis of the 21 cm hyperfine 
structure line of H I, where the WNM is detected in emission and the CNM in absorption. 
However, only upper limits of the temperature can be inferred because the contribution of 
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-3- 



nonthermal, or turbulent, broadening cannot be measured separately (Heiles 2001b). Ul- 
traviolet (UV) spectra of absorption lines can be used to make accurate temperature and 
turbulent velocity measurements of gas in the LISM by comparing the observed line widths of 
ions with different atomic masses. However, for distant lines of sight, the absorption features 
are typically severely blended, and most information regarding the influences of thermal and 
turbulent broadening mechanisms on the line widths is lost. Only crude estimates of the 
temperature can be made from the detection or nondetection of various ionization stages of 
observed elements (Welty et al. 2002). Therefore, nearby targets, and the intervening warm 
partially ionized gas, represent a unique opportunity to analyze unsaturated absorption 
lines with simple and well-characterized component structures so as to measure accurately 
the temperature and turbulence along lines of sight to these targets. 

After more than a decade of high resolution UV spectroscopy with the Hubble Space Tele- 
scope (HST), we are beginning to accumulate a critical number of observations of the LISM. 
Collections of observations presented by Redfield & Linsky (2002, 2004), and researchers ref- 
erenced therein, provide an opportunity to make fundamental physical measurements of the 
velocity, temperature, and turbulent velocity in the gas along the line of sight. Researchers 
investigating the LISM absorption lines towards nearby stars typically calculate a temper- 
ature and turbulent velocity for a given line of sight, but do not analyze the results in the 
context of the LISM as a whole. We present here the first systematic survey of temperature 
and turbulent velocity measurements for the LISM gas within 100 pc of the Sun. 



2. Calculating the Temperature and Turbulent Velocity 

By measuring the LISM absorption line widths of many ions, one can disentangle the 
contributions of thermal and turbulent broadening. The standard equation that relates the 
absorption line width, or Doppler parameter (6), and the temperature (T) and turbulent 
velocity (£) of the absorbing gas is 

b 2 = +£ 2 = 0.016629- + £ 2 , (1) 

m A 

where k is Boltzmann's constant, m is the mass of the ion, and A is the atomic weight of the 
ion in atomic mass units. Although only two line width measurements are required to solve 
Equation 1 for the temperature and turbulence, an accurate solution requires that the two 
ions have very different masses. Ideally, one would like to compare line widths of the two ions 
with the greatest mass differences that are observed in the LISM, hydrogen and iron. Iron 
is an excellent candidate because there are a number of resonance lines in the 2000-2600 A 
region that span a range of optical depths, permitting precise measurements of unsaturated 
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absorption lines. In addition, due to its large mass, thermal broadening is typically negligible 
and the line is intrinsically narrow, which improves the precision of the line fitting and 
facilitates the identification of individual velocity components. Redfield & Linsky (2002) 
presented an inventory of all measured Fe II LISM absorption line measurements toward 
stars within 100 pc. 

The light ions, however, pose more of a problem. Although the strong Ly-a line at 
1215 A is available for measuring the H I absorption line width, the line is heavily saturated 
and very broad, which makes estimating the stellar continuum difficult. This line can also 
be broadened by heliospheric and astrospheric absorption components (e.g., Linsky & Wood 
1996; Wood, Linsky, & Zank 2000). Deuterium is the next lightest ion, only twice as heavy 
as hydrogen. The D I absorption line is offset by roughly — 82kms~ 1 from the H I Lj-a line. 
For short lines of sight the D I line is typically unsaturated. Since the D I line usually lies 
on the extended wing of the broad stellar Ly-a emission line, it is relatively straightforward 
to interpolate the background continuum. D I is the subject of many papers and has been 
analyzed by many groups. Redfield & Linsky (2004) bring together the full collection of D I 
Ly-a observations towards stars within 100 pc with known velocity component structures, 
including many new measurements. 

In Figures la through lm, we present our measurements of the temperature and tur- 
bulent velocity for all 50 observed LISM components along 29 lines of sight to stars located 
within 100 pc that have absorption line width measurements of D I and an ion at least as 
heavy as Mg II. Sources for all of the Doppler width measurements used in Figures la-lm 
can be found in Redfield & Linsky (2002) and Redfield & Linsky (2004). For each veloc- 
ity component along a particular sightline, there are two plots to visualize the temperature 
and turbulent velocity measurement. On the left side is a color coded plot of temperature 
versus turbulent velocity. For a given measured Doppler parameter of a particular ion, Equa- 
tion 1 defines a curve in the temperature-turbulent velocity plane. The accompanying pair 
of dashed lines indicate the ± la error bars around the best fit curve (solid line). The black 
"x" symbol denotes the best fit value of the temperature and turbulent velocity based on 
the widths of the deuterium line and all other line widths measured from high resolution 
echelle (STIS E140H and E230H and GHRS Ech-A and Ech-B) spectra. The surrounding 
black contours enclose the la and 2a ranges for both values. The la and 2a error contours 
were determined by calculating the Ax appropriate for the 68.27% and 95.45% confidence 
level, respectively, as described in Bevington & Robinson (1992). Except for deuterium, all 
the other absorption lines are not fully resolved at moderate resolution. Therefore, system- 
atic errors, such as, inaccurate instrumental line spread functions, unresolved blends, and 
saturated lines, can potentially dominate the moderate resolution results. For this reason, 
the fits were made only to the deuterium measurement and all additional high resolution 



-5- 



measurements. The moderate resolution data are plotted as consistency checks. Two sight- 
lines, HR1099 and /3Cet, analyzed by Piskunov et al. (1997), use the H I Doppler width as 
opposed to D I. Although they did not provide an independent D I line width measurement, 
their H I Doppler width was forced to agree with the D I line. Therefore, despite the diffi- 
culty of analyzing the H I absorption line, these are likely to be reasonable estimates of the 
H I Doppler parameter. 

The plots on the right show another visualization for determining the temperature and 
turbulent velocity. Here, the observed Doppler parameter for each ion is plotted versus its 
atomic mass. The red and blue symbols represent the observations that were taken with 
high resolution and moderate resolution, respectively. The solid line indicates the best fit 
curve for the temperature and turbulence, as defined by Equation 1. The shaded region 
encompasses all fits within the la contour shown in the plot on the left. 

Both plots clearly demonstrate the need for using both the light ion, deuterium, and 
preferably the heaviest ion, iron, for calculating the temperature and turbulent velocity. 
The deuterium curve in the plots on the left is nearly vertical, whereas the curve for iron 
is practically horizontal, indicating that for deuterium thermal broadening dominates, and 
the line width is relatively insensitive to turbulent broadening. The reverse is true for iron. 
In the plots on the right, D I lies on the steep part of the curve, while practically all the 
remaining ions are on the flat part of the curve, where turbulent broadening dominates, 
and the line widths are independent of atomic mass and thermal broadening. Deuterium 
is significantly broader than all of the other ions and is the only serious constraint on the 
temperature of the gas. Even the next lightest ion observed in the LISM, C II, is already 6 
times the mass of deuterium, and predominantly horizontal in the plots on the left and on 
the flat part of the curve in the plots on the right. 

Although only two measurements are required to solve Equation 1, the intermediate 
mass ions provide important consistency checks against potentially unidentified systematic 
errors in the heavy and light ion measurements. Often, the lines of C II and O I have large 
errors in the Doppler parameter measurement because they are highly saturated. HZ 43 in 
Figure lj provides an excellent example of agreement among all of the measured ions. Only a 
handful of components show significant disagreement among the observed ions. In particular, 
Altair (Figures lc and Id) shows possible discrepancies in the observed C II Doppler width. 
This is most likely due to the C II absorption line being highly saturated. Low signal-to-noise 
observations or multiple closely spaced components can also make a precise Doppler width 
measurement difficult. Otherwise, this technique offers a unique opportunity to measure 
accurately the temperature and turbulent velocity in distinct collections of gas in the LISM. 

Table 1 lists the best fit temperature and turbulent velocity for all components shown 
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in Figures la through lm. The ± la error bars were determined by fixing each parameter 
in turn, and calculating the Ax appropriate for the 68.27% confidence level, as described in 
Bevington & Robinson (1992). To determine the best fit temperatures and turbulent veloci- 
ties for each velocity component, we weighted the observed Doppler parameters inversely to 
their measurement errors. We also include, for each component, the list of ions for which a 
Doppler width was included in the fit. Two is the minimum number of ions required for this 
technique, but our fits typically include Doppler widths from ~ 5 different ions, and in some 
cases as many as eight. Also included in Table 1, are the galactic coordinates and distance to 
the background star, and the weighted mean velocity measurements for the observed LISM 
absorption components. 



3. Temperature and Turbulent Velocity Distribution in the LISM 

The distributions of temperatures and turbulent velocities in the LISM are shown in 
Figure 2. The temperature distribution is strongly peaked, although significant outliers exist 
and range from K up to 12900 K. In particular, 18% of the sample, or between 8% and 38% if 
we include the la error bars, include measured cloud temperatures that are < 5000 K, which 
is in the thermally unstable regime (McKee & Ostriker 1977). Heiles (2001b) discussed a 
survey that is currently obtaining Zeeman splitting observations of H I 21 cm ISM absorption 
lines that may be able to provide a distribution of temperatures for more distant warm ISM 
environments. Our turbulent velocity distribution is likewise peaked, and the range of values 
extends from Okms^ 1 to 5.6kms _1 . The few very high turbulent velocity measurements are 
likely caused by unresolved blends. Because line widths can increase greatly when there 
are unresolved closely-spaced velocity components, which are common in the LISM (cf. 
Welty, Morton, & Hobbs 1996), undetected blends will increase the widths of all observed 
ions, although only by a small amount for deuterium. Therefore, the gas associated with 
a blended absorption feature will be interpreted as having an erroneously high turbulent 
velocity. 

Numerous circumstances specific to each observation can influence the precision of our 
temperature and turbulent velocity measurements. For example, a single absorption com- 
ponent along the line of sight, high resolution and high signal-to- noise observations, and a 
favorable H I line profile that leads to a straightforward continuum placement for D I, all 
contribute to a precise measurement of T and £. On the other hand, uncertainties in the 
derived values of T and £ are generally larger for lines of sight with several velocity compo- 
nents. Because of the range in precision of our measurements, we weight the mean inversely 
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to the measurement error, 



n 



Em/**) 

(T') = ^ , (2) 

E(iM 2 ) 

where Ti is the measured temperature for an absorbing cloud along a line of sight, Oi is the 
associated la error, and n is the total number of measurements, in this sample n = 50. In 
many cases the upper (erf') and lower (a^) bounds of the la error do not have the same 
magnitude. Therefore, when af ^ a~, we use the larger of the two quantities as <7j. The 
weighted mean temperature of the observed clouds in the LISM is ( T' ) = 6680 K. Likewise, 
the mean turbulent velocity weighted inversely by the measurement error can be calculated 
from Equation 2 by replacing T by £. The weighted mean turbulent velocity of clouds in the 
LISM is (O = 2.24 km s _1 . 

Although the temperature and turbulent velocity distributions shown in Figure 2 are 
significantly peaked about the weighted mean, the distributions are not consistent with 
a constant value of T or £ for all absorbers in the LISM. This is shown by calculating 
the xl goodness-of-fit measure under the assumption of a constant T and £ model for the 
LISM. The minimum value of x^(^lism) = 2.6 corresponds, as expected, to the weighted 
mean temperature, (T") = 6680 K, and likewise, the minimum value of x^(£lism) = 5.6 
corresponds to the weighted mean turbulent velocity, (£') = 2.24kms~ 1 . However, the 
large values of xl indicate that the observed variations of T and £ cannot be explained by 
the measurement errors alone, but instead the temperature and turbulent velocity must vary 
among collections of gas in the LISM. 

In order to quantify the observed variation in temperature and turbulent velocity in 
LISM absorbers, we calculate the weighted dispersion about the mean value, 
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as described by Bevington & Robinson (1992). Applying Equation 3 to our temperature 
measurements, we find that the dispersion about the weighted mean temperature is a it 1 ) — 
1490 K. By substituting £ for T in Equation 3, we find that the dispersion about the weighted 
mean turbulent velocity in LISM absorbers is cr/^\ = 1.03 kms -1 . Henceforth, any la error 
bars associated with temperature or turbulent velocity do not indicate the error in the 
weighted mean, as this assumes that the sample comes from the same parent population 
or equivalently that these quantities are constant in the LISM and our sample observations 
would be noisy measurements of that constant quantity. Instead, the la errors indicate 



the dispersion about the weighted mean. In Figure 2 the shaded distributions indicate 
subsets of the most precise measurements in our sample, those with the smallest measurement 
errors. Comparing the full distributions to the best measurement subset distributions clarifies 
the effect of weighting the mean and the dispersion about the mean by the inverse of the 
measurement errors. Since no information about the measurement errors is conveyed, the 
complete distribution is implicitly evaluated as if all measurements were made with the 
same precision. The best measurement subset distributions for both T and £ are more 
strongly peaked than the complete samples, because fewer outliers are included, thereby 
reducing the dispersion. Because the weighting favors the most precise measurements, the 
mean and dispersion about the mean closely correspond to the best measurement subsample 
distribution. 



3.1. Negative Correlation of Temperature and Turbulent Velocity 

In Figure 3 the temperature and turbulence measurements are compared in aggregate 
for all 50 absorbers. The subsample of filled symbols indicates a selection of the most precise 
measurements. The ± la error bars of some of our measurements can span almost the entire 
parameter space, and therefore do not provide much insight into the physical structure of the 
LISM along that particular line of sight. We have therefore selected a subsample of the most 
precise measurements, where the absolute error is less than the standard deviation of the 
entire sample, given above. This subsample includes more than a third of the entire sample. 
There appears to be a moderate negative correlation between the temperature and turbulent 
velocity. The correlation coefficient (r) is -0.47 for the entire sample and -0.35 for the precise 
subsample, with the probability that the distribution could be drawn from an uncorrelated 
parent population (P c ) is 0.010% and 3.5%, respectively (assuming no correlation of the 
errors). The probability, P c , increases when either \r\ or the sample size decreases. 

As is evident from the alignment of the error contours in Figure 1, the majority of 
temperature and turbulent velocity determinations show a correlation in their measurement 
errors. This covariance has approximately the same alignment as the negative correlation. In 
other words, irrespective of the presence or absence of a negative correlation, the covariance 
between the errors in T and £ will tend to smear out the measurements along a negatively 
sloped axis. To determine the significance of this effect on our sample, we calculate a large 
number of realizations (~ 10 5 ) for two different scenarios that have the same weighted mean 
and dispersion in temperature and turbulent velocity as our true dataset. A hypothetical 
dataset is used, which has the same number of elements as our true sample, but displays no 
correlation (r = 0). Each realization is a modification of the hypothetical dataset consistent 
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with the errors of our true sample. The two scenarios correspond to assumption that the 
errors are: (1) uncorrelated, where the axes of the error ellipses would be parallel to the T 
and £ axes, and (2) correlated, as is the case for our true sample, where the error ellipses are 
not aligned with the principal axes. The resulting distributions are shown in the bottom of 
Figure 3. The first scenario (solid curve) is a distribution of samples drawn from uncorrelated 
errors, which as expected, peaks at r ~ 0.0, and falls off dramatically at levels of moderate 
correlation (|r| ~ 0.3). The second scenario (dashed curve) also draws from our uncorrelated 
hypothetical sample but is varied according to the errors, where the errors are drawn from 
the true covariance matrix, taking into account the correlation of measurement errors. As 
expected, the distribution is slightly skewed to more negative correlations, peaking at r ~ 
—0.19, due to the smearing of the sample along a negatively sloped axis. The negatively 
sloping error ellipses in our sample lead to a shift of any distribution of correlations to slightly 
more negative values, thereby increasing the probability that the observed sample could be 
drawn from an uncorrelated parent population. In this case, the probability increases by 
about a factor of 10 for the second scenario compared to the first scenario where the errors 
are uncorrelated (see shaded regions in Figure 3), but the overall probability still remains 
relatively small (< 1%). Therefore, the covariance between temperature and turbulent 
velocity is a minor influence on the distribution of the sample, and the negative correlation 
between T and £ appears to be significant. 

A negative correlation between our temperature and turbulent velocity measurements 
implies that they may not be independent variables and that either some physical phe- 
nomenon, possibly pressure equilibrium, is causing the anticorrelation, or that systematic 
errors in our analysis are leading to the negative correlation. Possible sources of system- 
atic errors that are unique to each individual spectrum include the placement of the stellar 
continuum and the wavelength calibration. However, we expect these systematic errors to 
"randomize" within our sample, because they are different for each spectrum as a result of 
observing a variety of stars, where the location of the absorption is in a different part of the 
stellar continuum, and our use of different instruments with a range of spectral resolving 
powers. One important source of systematic error that could be ubiquitous in our sample, 
is the possible presence of unresolved blends. We attempt to minimize this problem by 
observing the narrowest LISM absorption lines at the highest available spectral resolution. 
However, one must consider the possible presence of several discrete absorbers with very 
similar projected velocities (cf. Welty, Morton, & Hobbs 1996) and their effect on our results 
when we assume that there is only one absorber along the line of sight. 

If unresolved blends are common in our sample, and multiple discrete absorbers are in- 
terpreted as single absorbers, then our line width measurements will be too broad. Because 
the unresolved blends will affect all resonance lines, independent of mass, the overestimate 
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of the line width for each ion leads to an overestimate of the turbulent velocity. Therefore, 
some of the high turbulence measurements are probably due to unresolved blends rather 
than to intrinsically high turbulent velocities. If we neglect the turbulent velocity measure- 
ments greater than 3kms _1 , for example, the negative correlation between temperature and 
turbulence is still present. The correlation coefficient is then r = —0.35 for the full sample 
and r = —0.44 for our best measurements subsample, and the probability (P c ) that the 
distribution results from an uncorrelated sample increases to 0.50% and 2.0%, respectively. 
The weak negative correlation may indicate a physical relationship between the temperature 
and turbulence in LISM clouds, perhaps involving pressure equilibrium. 



3.2. Temperature, Turbulent Velocity, and Cloud Mass 

The subtle connections between various cloud properties may help elucidate fundamen- 
tal characteristics of interstellar gas. In Figure 4, we compare the measured temperatures 
and turbulent velocities with a proxy for cloud mass, the D I column density. The D I 
measurements are taken from Redfield & Linsky (2004). As discussed in Section 2, a line 
width measurement of a light ion, in particular D I, is critical for an accurate determination 
of cloud temperature and turbulent velocity. Therefore, a D I column density is available for 
all of the sightlines discussed here. The deuterium-to-hydrogen (D/H) ratio appears to be 
constant within the Local Bubble, D/H= (1.56 ± 0.04) x 10~ 5 (Wood et al. 2004). There- 
fore, the deuterium column density can be used to estimate the hydrogen column density, 
which is itself a proxy for the mass of the cloud along the line of sight, if we assume that all 
these clouds have comparable densities. 

The plots in Figure 4 explore whether or not the temperature and turbulent velocity of 
the cloud vary with the cloud mass. Kinetic energy equipartition arguments suggest that, 
the turbulent motions of smaller clouds may be greater than those of more massive clouds. 
However, no correlation is discernible between column density and turbulent velocity, and 
equipartition does not seem to be important on the size scales of these clouds. We note that 
a single interstellar cloud may not have a unique turbulent velocity. Along a line of sight 
through an interstellar cloud, it is more likely that we are averaging over many turbulent 
environments within a single cloud structure, which could mask any equipartition signature. 
The presence of unresolved velocity structures along the line of sight could also disrupt 
the distribution of turbulent velocity measurements. No correlation is detected between 
temperature and column density either. Although the full sample seems to show a moderate 
positive correlation, that trend is not mirrored in the best measurement subsample, which 
shows no variation in T or £ with column density. These comparisons seem to indicate that 
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independent of cloud mass or size, the warm partially ionized clouds that populate the LISM 
share similar temperature and turbulent properties. 



3.3. Cloud Structure and the Routly-Spitzer Effect 

Implicit in our determination of the temperature and turbulent velocity for each compo- 
nent along a line of sight is the assumption that all the ions used in our analysis are members 
of the same collection of gas. Because all of the ions used in this paper are expected to be the 
dominant ionization stage, where the ionization fraction is > 50%, for physical conditions 
appropriate for the LISM (Slavin & Frisch 2002), this assumption is probably fair. Although 
some ions, such as, C II, Al II, Si II, and Fe II, clearly dominate the ionization distribution 
with > 95% of atoms expected in their ionization stage, other ions, including, D I, N I, O I, 
and Mg II, have less extreme ionization fractions ranging from 55%-85%. Wood et al. (2002) 
used LISM column density measurements along the line of sight toward Capella to make 
a direct comparison between observed ionization fractions of carbon, nitrogen, and silicon, 
and model predictions by Slavin & Frisch (2002). If differential ionization is acting such that 
various ions trace distinct interstellar medium, particularly since the degree of ionization 
fraction dominance spans a significant range for the observed ions, the measured line widths 
may not fit a single temperature and turbulent velocity. We can test the assumption that 
the observed ions reside in the same collection of gas by looking for systematic deviations in 
the Doppler width of each ion compared with the trend for other ions in the same component 
that would arise if individual ions were sampling different environments. For example, along 
distant lines of sight (~ lkpc), Routly & Spitzer (1952) noted that the velocity dispersion 
of the absorbers increased for clouds which had low Na I to Ca II ratios. This has been 
confirmed by more recent observations (Siluk & Silk 1974; Vallerga et al. 1993), and is at- 
tributed to large scale structures, such as old supernova shells, that can account for both 
the large velocity dispersion and the change in abundance ratio. It is possible that some- 
thing analogous to the Routly-Spitzer effect, whether due to varying depletion or ionization 
structure within an individual cloud, is acting on a much smaller scale, thereby separating 
the observed ions into slightly different environments within the cloud. We investigate this 
issue briefly here, while a detailed discussion of the ionization structure of the LISM will be 
provided in a subsequent paper. 

In Figure 5, we test for any systematic variations by showing the weighted mean de- 
viation for each ion from the best fit temperature and turbulent velocity solution for each 
cloud. The results are plotted against atomic mass as well as depletion for a typical diffuse 
cloud taken from the cool diffuse cloud along the ( Oph sightline (Savage & Sembach 1996). 
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Although the magnitudes of the depletions of the two diffuse clouds along the line of sight 
toward ( Oph are slightly offset from each other, the general trend in Figure 5 is the same 
regardless of which cloud is used as the depletion standard. The absorption features of both 
carbon and oxygen are typically saturated in our sample, which leads to their large dispersion 
about the weighted mean. Both deuterium and iron show very little deviation from the best 
fit value because they are at opposing ends of the range of atomic masses, and therefore, are 
the most influential in the best fit determination of the temperature and turbulent velocity. 
There is no indication of a systematic trend for different ions based on their depletion. In 
fact, the weighted mean of Ab from the best fit value is never greater than O.lkms -1 , which 
supports the assumption that these particular ions reside in the same collection of gas and 
provides accurate determinations of temperature and turbulent velocity. 



4. Estimate of the Thermal and Turbulent Pressures in the LISM 

From the distribution of the observed temperatures and turbulent velocities in the warm 
clouds in the LISM, we can estimate the mean thermal and turbulent pressures. The thermal 
pressure is defined as, 

P T = nkT, (4) 

where n is the number density of the particles in the gas, k is Boltzmann's constant, and T 
is the gas temperature. Likewise, the turbulent pressure is defined as, 

Pi = \p?, (5) 

where p is the mass density of the gas particles, and £ is the turbulent velocity of the 
gas. In order to calculate the gas density in the LISM, we assume cosmic abundances, or 
n(H)/n(He) ~ 10, and that the number of electrons equals the number of protons. The 
electron number density is assumed to be n c = O.ll^Qogcm -3 , as measured for the Capella 
line of sight through the Local Interstellar Cloud (LIC) using the C II* fine structure line 
(Wood & Linsky 1997), and using a nearby white dwarf (Holberg et al. 1999). The hydrogen 
number density for the LISM is assumed to be n(HI) = 0.1 cm" 3 (Redfield & Linsky 2000). 
The number density (n) is calculated from n = Yl rii, where % signifies the particles included 
in this computation: H I ions, electrons, protons, and helium ions. Similarly, the density (p) 
is calculated from p = Y n i m i- 

Given the mean temperature of the LISM, T = 6680 K and the measured dispersion of 
1490 K, the mean thermal pressure is calculated to be P^jk = 2280 K cm -3 with a dispersion 
about the mean of 520 K cm -3 , which is consistent with the value P/k = 1700 to 2300 cm -3 K 
estimated by Vallerga (1996) from Extreme Ultraviolet Explorer (EUVE) spectra. The mean 
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turbulent pressure and standard deviation, for a mean turbulent velocity of 2.24 km s" 1 and 
dispersion of 1.03 km s -1 , is Pg/k = 89 ± 82Kcm~ 3 . The thermal pressure thus dominates 
the turbulent pressure, Pt/P% ~ 26. Therefore, the small-scale motions in the warm clouds 
in the LISM, are dominated by thermal motions, as opposed to small-scale turbulence. The 
pressure estimates given above ignore the measurement error of n c in order to show the 
magnitude of the dispersion about the weighted mean pressures. If instead, we ignore the 
dispersion about the mean temperature and turbulent velocity, and include the measurement 
error associated with n c , the measured thermal pressure is P^/k = 2280^™ K cm -3 , and 
the turbulent pressure is P^/k = 89l 3 gKcm~ 3 . Different estimates of the hydrogen number 
density will not change this ratio significantly because the number density essentially cancels 
out of the ratio. 

The mean total pressure of the warm partially ionized clouds that populate the LISM 
is Ptot/k = Pr/k + P^/k = 2370 K cm -3 with a dispersion about he mean of 530 Kcm~ 3 . 
These clouds are surrounded by a larger structure, the Local Bubble, which presumably is 
characterized by hot (T ~ 10 6 K) low density (n c ~ 0.005 cm' 3 ) gas, which corresponds to 
a thermal pressure of P^/k ~ 11000 Kcm -3 (Sanders et al. 1977). The pressure imbalance 
between the warm and hot components of the LISM has been a persistent problem for 
understanding the structure of our local interstellar environment (Jenkins 2002). The present 
survey demonstrates that this discrepancy is common to all known local interstellar clouds 
within 100 parsecs, and that turbulent pressure is too small to make up the difference between 
the thermal pressure differences of the warm and hot LISM components. Recent results have 
stimulated a reexamination of the soft X-ray emission measurements of the Local Bubble, 
and hence the thermal pressure calculation of the hot medium. It turns out that 25-50% 
of the detected soft X-rays, initially wholly attributed to the hot gas in the Local Bubble, 
are the result of charge exchange between solar wind ions and interplanetary neutrals (Cox 
1998; Cravens 2000; Robertson & Cravens 2003). Therefore, the thermal pressure of the hot 
component of the LISM may be significantly lower than the initial estimates, and closer to 
agreement with the warm clouds of the LISM. 



5. Turbulent Mach Number 

We calculate the turbulent Mach number (M^) by taking the ratio of the observed 
turbulent velocity to the sound speed of the gas along the line of sight, 
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where the calculation of the number density (n) and mass density (p) are described in 
Section 4. The results of this calculation are shown as a function of cloud temperature in 
Figure 6. The gradual decrease in turbulent Mach number as a function of increasing cloud 
temperature that is seen in Figure 6 is a consequence of the increase in sound speed with 
temperature, as well as the anticorrelation of turbulent velocity with temperature, discussed 
in Section 3.1. Again, the best measurement subsample is displayed as red symbols, almost 
all of which are clustered at a turbulent Mach number of ~ 0.2. The weighted mean of the 
turbulent Mach number ((M/)) is calculated by replacing T by Mg in Equation 2, where 
(M^) = 0.19. Similarly, the weighted dispersion about the mean value (u(mm) is derived 
from Equation 3, where ct/m') — 0.11. The same data are also displayed in histogram form 
in Figure 7. Clearly, the macroscopic nonthermal motions of warm clouds in the LISM are 
significantly subsonic, in contrast to the highly supersonic turbulent motions in cold neutral 
clouds (Heiles & Troland 2003), and in agreement with optical observations of cool clouds 
(Welty, Morton, & Hobbs 1996; Dunkin & Crawford 1999). 



6. Spatial Distribution of Temperature and Turbulent Velocity 

In Figures 8 and 9, the observed temperature and turbulent velocity, respectively, for 
each line of sight are shown in Galactic coordinates. The size of a symbol is inversely pro- 
portional to the target's distance, and the color of a symbol indicates the observed quantity, 
as shown by the scale at the bottom of each plot. Sightlines that contain more than one 
absorber, or sightlines to multiple targets, such as binary systems, include measurements of 
all clouds by separating the shadings with a dotted line. In this way, each sightline will have 
between one and three different shadings to signify up to three absorbers along the line of 
sight. 

In a number of cases, both the temperature and turbulent velocity seem to be spatially 
correlated. For example, the closest pair of sightlines, a Cen A and a Cen B, have almost 
identical temperature and turbulent velocity measurements. The same is true for another 
close pair, Capella and the more distant white dwarf G191-B2B, which are only separated by 
7°. The +20kms~ 1 component is common in the spectra of both stars, and the individual 
temperature and turbulent velocity measurements are practically identical. On the other 
hand, a closely grouped trio of stars in the direction of the North Galactic Pole, share only 
moderately similar cloud characteristics. With angular separations of 5°-8°, the observed 
spectra of 31 Com, HZ 43, and GD 153, each show one absorption component with roughly 
the same velocity. The derived temperature measurements are similar, ranging from 7000- 
8200 K, but the turbulent velocity measurements, while all less than the LISM mean, range 
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from no detected turbulent velocity to 1.7kms _1 but all are formally consistent with £ = 
O.Okms -1 . The temperatures and turbulent velocities for the two closely spaced velocity 
components in the lines of sight to Procyon and 61 Cyg A are identical because the line 
widths were assumed to be the same due to the blending of the absorption components. 

In order to quantify the magnitude of the spatial correlations in our measurements, we 
have compared the absolute values of the temperature and turbulent velocity differences for 
all unique pairings of the stars listed in Table 1. The results are plotted in Figure 10 as 
a function of the angular separation of the pair of sightlines, in 10° bins. We also applied 
this procedure to the absolute values of the observed projected velocity difference (Av) of 
the absorbing cloud and the total pressure, which is a function of both temperature and 
turbulent velocity. 

It has been known for some time that the projected velocity of absorbing material in the 
LISM shows a strong spatial correlation (Crutcher 1982). In fact, single three dimensional 
velocity vectors have been used to describe the bulk flow of the closest interstellar clouds, 
and is one of the predominant reasons that LISM absorbers are thought of as roughly homo- 
geneous cloud-like entities, rather than chaotic and filamentary objects (Lallement & Bertin 
1992; Lallement et al. 1995). The spatial correlation in velocity is clearly seen in the top 
panel of Figure 10, where both the absolute velocity difference, and the dispersion about the 
weighted mean, are significantly less for sightlines with angular separations < 20°. The dot- 
ted lines indicate the weighted mean of measurements from 0°-20°, and then from 20°-60°. 
The absolute value of the weighted mean velocity separation for close pairs is 21% that of 
more distant pairs. The same trend seems to be true for temperature, where the absolute 
value of the weighted mean temperature difference of pairs with separations < 20° is 57% of 
the weighted mean of more distant pairs. For turbulent velocity measurements the difference 
between the closest and furthest pairings is only marginal, just 81%. The total pressure of 
the cloud shows a significant spatial correlation, where the weighted mean of pairings < 20° 
apart is 50% of the weighted mean of pairings with greater angular separations. The spatial 
correlation of total pressure may be related to the negative correlation between tempera- 
ture and turbulent velocity discussed in Section 3.1. Although temperature and turbulent 
velocity may each vary over short angular distances, pressure equilibrium may be acting to 
minimize the total pressure variation over small angular distances. 

Temperature should be useful as a unique cloud characteristic, because like velocity the 
total measurement range is large but over small angular separations the mean temperature 
range is relatively small. However, the correlation of temperature with small angular separa- 
tions is not as tight as for velocity. The turbulent velocity measurement may have less to do 
with the gross characteristics of a specific interstellar cloud, but could vary greatly inside a 
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cloud from stratification of the particular physical structure of the gas due to external forces 
such as cloud-cloud interactions. 



7. Conclusions 

We present here the first systematic survey of temperature and turbulent velocity in 
the warm partially ionized clouds of the LISM. The collections of observations presented by 
Redfield & Linsky (2002, 2004), and researchers referenced therein, provide an opportunity 
to make fundamental physical measurements of the LISM, including 50 independent tem- 
perature and turbulent velocity measurements of the gas along 29 lines of sight in the LISM 
within 100 pc. The results of this work can be summarized as follows: 

1. We find that in order to make an accurate and precise measurement of the gas tem- 
perature and turbulent velocity, the line widths of deuterium, and an ion at least as 
heavy as magnesium are required. Without deuterium there is no serious constraint on 
the temperature, and without a relatively heavy ion, there is no serious constraint on 
the turbulent velocity. High resolution observations are absolutely essential, because 
except for deuterium, the absorption lines are not resolved at moderate resolution, and 
systematic errors can dominate the measurements, leading to erroneous results. 

2. The distributions of temperature and turbulent velocity of LISM gas located are both 
significantly peaked. The weighted mean and standard deviation for the temperature 
is 6680 ± 1490 K. The weighted mean and standard deviation for the turbulent velocity 
distribution is 2.24 ± 1.03 km s -1 . The high turbulent velocity measurements are most 
likely due to unresolved blends. The standard deviation indicates the range of variation 
in the LISM, rather than the measurement error of the mean value. We find a significant 
number (18%, or between 8% and 38% including the la errors) of sightlines with 
temperatures < 5000 K, below the stable WNM temperature regime defined by McKee 
& Ostriker (1977). 

3. A moderately negative correlation is detected between temperature and turbulent ve- 
locity. The correlation coefficient (r) is -0.47 for the entire sample and -0.35 for the 
precise subsample, with the probability that the distribution could be from an uncor- 
rected parent population (P c ) of 0.010% and 3.5%, respectively. The covariance of the 
errors of the temperatures and turbulent velocities do not dominate the determination 
of this negative correlation. Pressure equilibrium among the warm clouds may be the 
source of this anticorrelation. 
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4. The small-scale thermal motions dominate the turbulent motions. For the mean tem- 
perature and turbulent velocity, we calculate a mean thermal pressure of P^/k = 
2280 Kcm -3 and a dispersion of 520 K cm -3 , which is consistent with the value P/k = 
1700 to 2300 cm -3 K estimated by Vallerga (1996) from EUVE spectra. The mean ther- 
mal pressure dominates the mean turbulent pressure, P^/k = 89 K cm -3 with a disper- 
sion about the mean of 82Kcm -3 , by more than an order of magnitude, Pt/P% ~ 26. 
Turbulent pressure in LISM clouds cannot make up the difference in the apparent pres- 
sure imbalance between warm LISM clouds and the surrounding hot gas of the Local 
Bubble. 

5. We find that the turbulent motions of the warm partially ionized clouds in the LISM 
are significantly subsonic. The weighted mean and standard deviation of the turbulent 
Mach number for the entire sample is 0.19 ± 0.11. Although some turbulent Mach num- 
bers approach the sound speed, due to relatively low temperatures of the clouds, these 
temperatures are not accurately known. The best measurement subsample clusters 
almost exclusively around turbulent Mach numbers of ~ 0.2. 

6. We have calculated the spatial distribution of the temperature and turbulent velocity 
for the 50 velocity components in the LISM. These individual measurements can now 
be combined with additional physical properties of LISM clouds, such as projected 
velocity and depletions, to determine the three-dimensional morphological structure of 
gas within 100 pc. This work will be presented in subsequent papers. 
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Fig. la. — Determination of the temperature (T) and turbulent velocity (£) for each LISM 
absorption component observed within 100 pc. These measurements require an observation 
of a light ion such as D I, and at least one other ion at least as heavy as Mg II. Each LISM 
component is labeled in the plot on the left with the sightline designated by the background 
star, and the velocity of the LISM absorption. The plot on the left shows the best fit 
Doppler parameter (solid lines) and ± la error bars (dashed lines) for each ion observed, 
colored coded and labeled toward the right side of the plot. Each curve is determined by the 
measured Doppler parameter and mass of the observed ion, as described in Equation 1. The 
"x" symbol indicates the best fit temperature and turbulent velocity given the widths of all 
observed ions, while the surrounding black contours indicate the ± la and ± 2<r error bars. 
The plot on the right shows the measured Doppler parameters as a function of the atomic 
mass of the observed ion, which are specifically indicated in the top axis. Red symbols 
indicate that the LISM absorption line was measured from a high resolution observation and 
used to make the fit, and blue symbols indicate that the measurements were from moderate 
resolution observations and are used as a consistency check. The best fit temperature and 
turbulent velocity is given in the right plot, and the best fit curve is displayed by the solid 
line, and the shading includes the ± la range of possible outcomes given the derived T and 
£ parameters and their respective errors. 
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Fig. 2. — Distributions of measured temperatures and turbulent velocities in the LISM 
for clouds within 100 pc. The temperature distribution is peaked at the canonical LISM 
temperature ~ 7000 K. The weighted mean of our sample is ( T" ) = 6680 K, while the dis- 
persion about the mean is 1490 K. The weighted mean of the turbulent velocity distribution 
is (O =2.24 km s -1 , with a dispersion about the mean of 1.03 km s -1 . The weighted mean 
is indicated by a dotted line. The shaded subsamples represent our most precise measure- 
ments, defined by their absolute error being less than the standard deviation of the entire 
sample. 
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Fig. 3. — The top plot shows the turbulent velocity as a function of temperature for our 
complete sample of LISM absorbers within 100 pc. The filled symbols indicate our most 
precise measurements, defined by their absolute error being less than the standard deviation 
of the complete sample (see Figure 2). The dashed lines indicate the weighted mean values 
of the temperature and turbulent velocity. The bottom plot shows the distribution of cor- 
relation coefficients in ~ 10 5 realizations of two different scenarios. The solid curve shows 
the distribution of a hypothetical dataset that shows no correlation, varied over the true 
errors, but assuming the errors in temperature and turbulence are independent and uncorre- 
lated. The dashed curve indicates the same hypothetical sample varied according to the true 
error bars, but now taking into account the actual alignment of the error contours, which 
indicate a correlation in the errors of the variables (see error ellipses in Figure 1). For all 
realizations, the weighted means and dispersions of both variables (T and £) were identical 
to the true sample. The vertical line indicates the measured correlation of our true sample. 
Although the probability that this value of correlation is a result of an uncorrelated parent 
population increases when we take into account the correlated errors (see shaded regions), 
the overall probability is still very low (< 1 %). Therefore, the covariance between the errors 
of temperature and turbulence do not appear to dominate the distribution of the sample, 
and the moderately negative correlation between temperature and turbulence seems to be 
significant. 
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Fig. 4. — Cloud temperature and turbulent velocity are plotted against D I column density. 
The filled symbols indicate the best measurement subsample. The dashed lines indicate 
the weighted mean values of the temperature and turbulent velocity. The D I column 
density can be converted to a H I column density by applying the locally constant ratio 
D/H= (1.56 ± 0.04) x 10 -5 (Wood et al. 2004). The hydrogen column density can be used 
as a proxy for the cloud mass. 
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Fig. 5. — Departures of the weighted mean values of the Doppler parameters are compared 
with the best fit values for each individual ion. The top plot shows the weighted mean devi- 
ation against atomic mass, and the bottom plot shows the same against a typical depletion 
for each ion in the diffuse interstellar medium, taken from the cool diffuse cloud along the 
C Oph sightline (Savage & Sembach 1996). Although the magnitudes of the depletions of 
the two diffuse clouds along the line of sight toward ( Oph are slightly offset from each 
other, the general trend displayed here is the same regardless of which cloud is used as the 
depletion standard.. Carbon and oxygen have large deviations and uncertainties about the 
weighted mean because their absorption features are typically saturated, which inhibits a 
precise measurement of the Doppler parameter (b). No systematic trends are evident for this 
class of ions, supporting the assumption that they are all tracing the same collection of gas 
in each velocity component. 
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Fig. 6. — A plot of the turbulent Mach number as a function of cloud temperature. The 
gradual decrease in turbulent Mach number as a function of increasing cloud temperature is 
a consequence of the increase in sound speed with temperature, as well as the anticorrelation 
of turbulent velocity with temperature, discussed in Section 3.1. However, the vast majority 
of measurements, in particular the best measurement subsample, indicate a turbulent Mach 
number ~ 0.2, shown by the dashed line. A Mach number of unity is indicated by the dotted 
line. 
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Fig. 7. — Distribution of the observed turbulent Mach numbers in the LISM. The weighted 
mean of our sample is (Ml) =0.19, while the dispersion about the mean is 0.11. The 
weighted mean is indicated by a dotted line. The best measurement subsample is indicated 
by the shaded distribution. The high values of the turbulent Mach number may be spurious 
due to the possible presence of unresolved closely spaced velocity components. Turbulence 
in clouds within 100 pc is clearly subsonic. 
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Fig. 8. — The spatial distribution of LISM cloud temperatures is shown in Galactic coor- 
dinates. The size of each symbol is inversely proportional to the distance, and the color of 
the symbol indicates the temperature of a cloud along the line of sight, as given by the scale 
at the bottom of the plot. If more than one absorbing cloud lies along the line of sight, 
additional cloud temperatures are indicated by nested circles separated by a dotted line. 
This technique is also used for measurements of different members of a single multiple star 
system, such a Cen. The temperatures for the two velocity components along the lines of 
sight to Procyon and 61 Cyg A are identical because the line widths of the closely spaced 
velocity components were assumed to be the same. 
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Fig. 9. — Similar to Figure 8, but displaying the turbulent velocity of LISM material within 
100 pc. 
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Fig. 10. — Differences between the measured physical parameters, velocity, temperature, 
turbulent velocity, and total pressure, as a function of angular separation. The bin size is 
10°, and the weighted mean for each bin is depicted by the circle symbol, while the dispersion 
about the weighted mean is indicated by the error bars. On average, 20 pairings are used 
in each bin, while the minimum, in the 10°-20° bin, is 6 pairings. A spatial correlation 
is obvious for the velocity measurements, where the absolute velocity differences, and the 
magnitude of the dispersion about the mean are much smaller for pairs of sightlines that are 
separated by 20° or less. The weighted mean for separations from 0°-20° and 20°-60° are 
indicated in each plot by the dotted line. A similar spatial correlation seems to be evident 
in the total pressure plot, and to a lesser degree in the individual temperature and turbulent 
velocity plots. 
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Table 1. Individual Temperature and Turbulent Velocity Measurements 



III) 


Other 


; 


b 


# 


Name 


(°) 


(°) 


128620 


a Ccn A 


315.7 


-0.7 


128621 


a Con B 


315.7 


-0.7 


22049 


e Eri 


195.8 


-48.1 


201091 


61 Cyg A 


82.3 


-5.8 



distance 
(pc) 



{v) 
(km s _1 ) 



T 
(K) 



(km s _1 ) 



26965 
165341 



197481 
62509 



Procyon 



40 Eri A 
70 Oph 



200.8 -38.1 
29.9 11.4 



155886 36 Oph A 358.3 6.9 

131156A £ Boo A 23.1 61.4 

39587 x 1 Ori 188.5 -2.7 

20630 k 1 Cet 178.2 -43.1 



AU Mic 
Gem 



C Dor 



12.7 
192.2 



-36.8 
23.4 



v Peg 



Cap 



G191-B2B 156.0 



GD 153 317.3 84.8 

111812 31 Com 115.0 89.6 



1.3 
1.3 
3.2 
3.5 



5.0 
5.1 



5.5 
6.7 
8.7 

9.2 



9.9 

10.3 



34029 


Capclla 


162.6 


4.6 


12.9 


432 


Cas 


117.5 


-3.3 


16.7 


11443 


a Tri 


138.6 


-31.4 


19.7 


22468 


HP, 1099 


184.9 


-41.6 


29.0 


4128 


Cet 


111.3 


-80.7 


29.4 


120315 


Alcaid 


100.7 


65.3 


30.9 




HZ 43 


54.1 


84.2 


32.0 


82210 


DK UMa 


142.6 


38.9 


32.4 


62044 


a Gem 


191.2 


23.3 


37.5 



70.5 
94.2 



-18.45 

-18.14 

18.73 

-3.0 

-9.0 

23.0 

19.76 

21.73 

-26.50 

-32.53 

-43.34 

-17.1 

-20.9 

-25.02 

-28.40 

-17.69 

23.08 

20.84 

13.35 

7.36 

-21.45 

31.84 

19.65 

13.9 

8.41 

21.48 

9.15 

17.89 

13.65 

21.90 

14.80 

8.20 

9.14 

1.63 

2.6 

-3.02 

-6.52 

9.41 

32.26 

21.77 

8.8 

1.73 

-7.48 

-2.22 

-12.06 

-20.48 

19.19 

8.61 

-5.04 

-3.37 



± 0.32 

± 0.08 
± 0.65 
± 1.4 
± 0.7 
± 1.0 
± 0.75 
± 0.22 
± 0.06 
± 0.32 
± 0.14 
± 1.2 
± 1.2 
± 0.84 
± 0.58 
± 0.60 
± 0.66 
± 0.20 
± 0.15 
± 0.46 
± 0.34 
± 0.60 
± 0.76 
± 1.3 
± 0.82 
± 0.51 
± 0.68 
± 0.57 
± 0.75 
± 0.04 
± 0.02 
± 0.02 
± 0.21 
± 0.22 
± 3.4 
± 0.51 
± 0.39 
± 0.61 
± 0.15 
± 0.16 
± 1.2 
± 0.39 
± 0.43 
± 0.19 
± 0.49 
± 0.45 
± 0.15 
± 0.38 
± 0.17 
± 0.37 



+ 1200 

- 1100 
+ 330 

- 320 
-I- 8(H) 

- 830 
± 880 
± 880 
+ 660 

- 630 
■I 000 

- 630 
± 450 
+ 3000 

- 2300 
+ 2100 

- 1700 
± 2100 
+ 2000 

- 2200 
± 2400 
+ 2700 

- 2400 
± 560 

± 830 
+ 730 

- 680 
+ 1900 

- 1700 
I 2900 

- 2200 
± 2700 
± 1200 
+ 3100 

- 2600 
+ 1600 

- 1500 
+- 3500 

- 3000 
+ 2300 

- 2100 
+ 1400 

- 1300 
+ 800 

- 880 
+ 3100 

- 2600 
+ 3900 

- 3400 
± 1500 
+ 800 

- 1100 
± 1400 

6300 ± 2900 
12400 ± 2800 
n +4400 
u - 

sarin + 2500 
8900 _ 2300 

7500 + 2100 
6750 ± 240 

7200 + 100 ° 
/zuu _ 120Q 

8600 ± 1600 

1000 + 1900 
iuuu _ 1000 

1700 +110 ° 

ofinn + 4400 
3600 Z 3000 

5500 + 10600 
oouu - 5500 

1 2900 + 3800 
I2yuu _ 3300 

11700 + 4100 
iuuu _ 3600 



5100 
5500 
7410 
6850 
6850 
6710 
6710 
8120 
2700 
1700 
3300 
12300 
12600 
12500 
5870 
5310 
7000 
5200 
3600 
5800 
8700 
6100 
9000 
7000 
7700 
6700 
9760 
7700 
8900 
7900 
8800 
7100 



0200 



+ 1400 
- 1300 

4400 +280 ° 
Muu - 2400 

7000 +290 ° 
/uuu _ 2800 

soon + 100 ° 

8200 _ 1400 



121 + 033 

i ' zl - 0.49 

-, ,, 7 +0.34 

1 - i ' - 0.41 

2 03 +°' 41 
/ ' UJ - 0.45 

2.08 ± 0.64 
2.08 ± 0.64 
-i 01 + 0-35 
'-■^ 1 -0.45 
1 01 + 0.35 
i '^ i - 0.45 

5 + 1 - 2 
"■•> -0.5 

3 64 +°' 42 
,3,D * - 0.44 

3.3 ± 1.1 

2.31 ± 0.37 
+ 1.2 

- 0.0 
+ 0.90 

- 0.63 
+ 0.7 

- 1.4 
+ 0.46 

- 0.51 
1.68 ± 0.23 

+ 0.15 

- 0.17 
+ 0.28 

- 0.32 
, 17 +0.34 
z1 ' - 0.52 
1.48 ± 0.92 
4.30 ± 0.93 

+ 0.79 

- 0.59 
+ 0.27 

- 0.32 

5 47 + 039 
°'*' -0.41 
+ 0.38 

- 0.48 
+ 0.32 

- 0.39 
+ 1.1 

- 0.0 
+ 1.7 

- 0.0 
+ 1.7 

- 1.3 

1.18 ± 0.47 

n nn + °- 85 
u.uu _ 00 

2.30 ± 0.25 

1.31 ± 0.76 

2.29 ± 0.44 

5.6 fA 

134 + °' 24 

1 '°* - 0.31 

17 +0.8 

1 -' - 1.7 

1 35 + 018 

00 + 1 - 1 

uu -0.0 
2.46 ± 0.45 

3 46 + 0131 
C1 '* D - 0.63 

3.93 ± 0.22 

+ 0.7 



0.63 

1.4 

2.33 



2.38 
2.64 



1.93 

1.67 



2.34 
1.68 

1.3 



1.7 

3.7 

1.58 

3.82 

1.78 

3.27 

1.2 

0.00 



- 1.5 
+ 0.7 

- 1.1 
+ 0.56 
- 0.89 

+ 0.37 

- 0.44 
+ 0.40 

- 0.51 
+ 0.37 

- 0.39 
+ 1.1 

- 1.2 
+ 0.99 

- 0.00 



CII, NI, OI, MgH, A1II, Sill, Fell 

Mgll, Fell 

CII, OI, Mgll, Sill, Fell 

Mgll 

Mgll 

Mgll, Fell 

Mgll, Fell 

Mgll 

CII, OI, Mgll, Sill, Fell 

CII, OI, Mgll, Sill 

CII, OI, Mgll, Sill 

CII, Mgll, Fell 

CII, Mgll, Fell 

CII, Mgll, Fell 

Mgll, Fell 

CII, OI, Mgll, Sill 

CII, Mgll, Sill, Fell 

CII, OI, Mgll, Sill, Fell 

CII, OI, Mgll, Sill, Fell 

CII, OI, Mgll, Sill 

Mgll 

OI, Mgll, Fell 

OI, Mgll, Fell 

CII, OI, Mgll, A1II, Sill, Fell 

CII, OI, Mgll, A1II, Sill, Fell 

CII, NI, OI, Mgll, A1II, Sill, Fell 

Mgll, Fell 

Mgll, Fell 

Mgll, Fell 

CII, OI, Mgll 

CII, OI, Mgll 

CII, OI, Mgll 

CII, NI, OI, Mgll, A1II, Sill 

CII, NI, OI, Mgll, A1II, Sill 

CII, Mgll, Sill 

CII, OI, Mgll, Sill, Fell 

CII, NI, OI, Mgll, Sill, Fell 

CII, OI, Mgll, A1II, Fell 

Mgll 

Mgll 

CII, OI, Mgll, Sill, Fell 

CII, OI, Mgll, A1II, Sill, Fell 

CII, OI, Mgll, Sill, Fell 



CII, NI, OI 
CII, NI, OI 
CII, NI, OI 
CII, NI, OI 
CII, NI, OI 
CII, NI. OI 



Mgll, A1II, Sill, Fell 
Mgll, A1II, Sill, Fell 
Mgll, Sill, Fell 
Mgll, A1II, Sill, Fell 
Mgll, A1II, Sill, Fell 
Sill, Fell 



CII, OI, Mgll, Fell 



Note. — This table includes only those measurements that include independent fits to the line width of deuterium and an ion at least as heavy as 
magnesium. 



